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Abstract 

In this paper, we present a numerical scheme to solve the initial-boundary value 
problem for backward stochastic partial differential equations of parabolic type. Based 
on the Galerkin method, we approximate the original equation by a family of backward 
stochastic differential equations (BSDEs, for short), and then solve these BSDEs by the 
time discretization. Combining the truncation with respect to the spatial variable and 
the backward Euler method on time variable, we obtain the global L 2 error estimate. 
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1 Introduction 

Let T G (0, -|-oo), (12, Y, F, P) be a complete probability space and F = { JJ, t G [0, T]} 
be the natural filtration generalized by a 1-dimensional Wiener process {W(t) : t G [0, T]} 
satisfying the usual conditions. The purpose of this work is to present a numerical scheme 
for solving the following backward stochastic parabolic differential equation (BSPDE, for 
short): 

{ dq(t, x ) = (—A q(t, x ) + f(t, x, q(t , x),r(t, x)))dt + r(t, x)dW (£), in [0, T) x D, 

q(t,x ) = 0, on [0, T) x dD , (1.1) 

q(T,x) = qr(x), in D. 
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Here, D C R is a bounded domain with C 2 boundary (d e N), A = Ylk=i^ 2 /^ x l i s the 
Laplacian, and / and qr are given data satisfying suitable conditions to be given later. 

BSPDEs are nontrivial extensions of BSDEs, which possess interesting theoretical val¬ 
ues and were originally introduced in the study of optimal control problems as the adjoint 
equation appeared in the Pontryagin maximum principle when the controlled system is a 
stochastic parabolic differential equation (see, e.g., [11, 15]), then are proved usefully in 
nonlinear filtering [20], mathematical finance [5] and so on. 

The well-posedness of general BSPDEs driven by Wiener process has been considered in 
a number of papers using either duality techniques [15, 30] or martingale representation and 
fixed point arguments [12, 22], 

In recent years, the study of numerical solutions to stochastic differential equations be¬ 
comes an active topic, and has attracted considerable attention in many fields such as control 
theory and mathematical finance. 

Up to now, some numerical schemes for the forward stochastic partial differential equa¬ 
tions have been presented: the Galerkin approximation ([7]), the finite difference method in 
space and time (see, e.g., [8, 21]), the finite element method (see, e.g., [23, 26]), the stochastic 
Taylor expansion method (see, e.g., [13]), the Wiener chaos expansion method (see, e.g., [9]) 
and so on. 

On the other hand, there are also several algorithms for solving BSDEs. First of all, based 
on the four step scheme [17], and using the relations between BSDEs and PDEs, Douglas 
et al. [4], Milstein and Tretyakov [19] obtained a numerical algorithm to solve BSDEs. The 
second kind of algorithms is the backward Euler method (see [1, 3, 6, 10, 29] for more details), 
for which the most difficult part is to calculate the conditional expectation. The third one 
is the random walk approach, in which the Brownian motion is replaced by a scaled random 
walk (see [16] and the references therein). We should also mention works by Wang and Zhang 
[24] presenting the finite transposition method, Lototsky et al. [14], Briand and Labart [2] 
giving the Wiener chaos expansion method and so on. 

Compared with the development of numerical methods for BSDEs and forward stochastic 
partial differential equations, the study of numerical schemes for BSPDEs is quite limited. To 
the author’s best knowledge, there exists no published work in this direction. Here we should 
mention the work of Yannacopoulos et al. [27], which only listed an idea on the numerical 
scheme for solving BSPDEs. Their method depends on the Wiener chaos expansion, however 
they did not prove the convergence speed. 

The aim of this work is to provide a numerical scheme for solving Eq. (1.1) based on 
the semidiscrete Galerkin approximation. This scheme is divided into two steps. Firstly, we 
approximate Eq. (1.1) by a family of BSDEs; then, adopting the backward Euler method, 
we give numerical solutions to the related BSDEs. This work is an improved and complete 
version of [25, Chapter 5]. 

Compared to the Wiener chaos expansion method listed in [27], it seems that the Galerkin 
scheme is easier to operate, at least for some special cases. Indeed, if the domain is “good” 
enough and the eigenvalues, eigenfunctions are easy to be computed, we can take the knite- 
dimensional approximation spaces to be the ones spanned by suitable eigenfunctions. 

The rest of the paper is organized as follows: In Section 2, we introduce our general setting 
and review the well-posedness of Eq. (1.1). In Section 3, we make use of the space-discretised 
Galerkin approximation to Eq. (1.1) and construct the desired finite-dimensional spaces; 
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then Eq. (1.1) is approximated by a family of BSDEs. We prove the strong convergence 
in appropriate spaces and obtain the rate of the convergence in the L 2 norm with respect 
to the space variable. In Section 4, we adopt the backward Enler method, for which the 
convergence and error analysis are also provided. 

2 Preliminaries 

Let H be a Hilbert space with norm || • ||# and inner product (•, • ) H . The following classes 
of processes will be frequently used throughout this paper: 

• For any t G [0, T\, H) is the space of all JF r measurable* //-valued random 

variables £ satisfying ||£|| 2 2 = EIISIIh < °o; 

• Cf([ 0, T]; L 2 (f2; //)) is the space of all F-adapted, continuous, //-valued stochastic pro¬ 
cesses X satisfying H^H^ao.T];^^)) = su Pte[o,T] e ITO)IIh < oo; <^([0, T];//)) 

is the subspace of CVQO, T]; L 2 (fl; //)) such that E(sup tg [ 0;T ] ||A"(/)||^) < oo; 

• Lp(0 x (0, T); H ) is the space of all F-adapted, //-valued stochastic processes Y satis¬ 
fying ||y||^ 2 (nx(0)T);H) = E(/ 0 T \\Y(t)\\Ut) < oo; 

• D k,p (H) (k , p > 1) is the space of all TV-rneasurable, H® k -v, alued, A;-times Malliavin 
differentiable random variables £ satisfying 

k 

ncin.p := (E(i€i p )+E E (ii^ir H «)) 1/P < +°°; 

3 =i 

• M 2,p (p > 2) is the space of all TV-measurable square integrable random variables £ 
which admits a stochastic integral representation: 

£ = E£ + f T u(t)dW (t), 

Jo 

where u(-) is a progresively measurable process satisfying sup 0<t<r E|u(/)| p < +oo; 

• L„ 2 (//) is the space of all //-valued progressively measurable processes {w(/)}te[o,r] 
satisfying 

(i) For almost all t G [0,T], u{t) G 0 1,2 (/Z); 

(ii) E(/ q T \u(t)\ 2 dt + / q T / q T \Dgu(t)\ 2 d,9dt) < oo. 

Besides, denote by (•, •) and | • | respectively the inner product and norm in different Eu¬ 
clidean spaces, which can be identified from the context; by A* the transport matrix of 
A. 

Let us introduce the following two assumptions: 

(Al) / : R + x D xR d xR d —>■ R d is Holder continuous with respect to t, i.e., 
there exists a positive constant L, such that \f(t\,x,y,z) — f(t 2 ,x,y,z)\ < Ly/\t i — ^21 7 f or 
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any ti, t 2 G R + , x G D, y, z G H d ; and has continuous and uniformly bounded first and 
second partial derivatives with respect to y and z (still denote this bound by L). Moreover, 
/(■>■> 0, 0) G Z/ 2 (0, T; Hq(D)). 

(A2) q T G Hq(D)) n B 2 *(L 2 (D)), 

E| DoC[t — Do'Qt\ ^ iW|0 — O '|, 

max |E|g T |^ 2(D) , sup El-Deg^l 9 , sup sup E|L>^R e g T | 9 ) < M, 
ee[o,T] 0'e[o,T] 0e[o,T] ' 

where q > 4 and M is a positive constant. 

Define A : D(A) = H 2 (D) C\Hq(D) —» L 2 (Z1) by Af = A/, for all / G D(A). It is easy 
to show that A is the infinitesimal generator of a Co-semigroup {e At }t> o on L 2 (D). Then 
Eq. (1.1) can be re-writen as the following abstract form: 


( 2 . 1 ) 


f dq{t) = ( - Aq(t) + f{t, -> f q{t),r{t)))dt + r(t)dW(t), t G [0 ,T), 
\ q(T) = q T - 


Now, we recall the definition of solution to Eq. (2.1) and its well-posedness. We refer 
the reader to [12, 22] for details. 

Definition 2.1. A pair of random fields ( q,r) G Cp([0, T]; L 2 (D; L 2 (D))) xL|(Dx(0, T); L 2 (D)) 
is called a mild solution to Eq. (2.1), if for every t G [0, T\, it holds that 

r T r T 

q{t)+ / e A ^f(s,-,q(s),r(s))ds + / e A ^r(s)dW(s) = e^V, a.s. 

Jt Jt 

Theorem 2.1. Suppose that (Al) holds and the terminal condition qT G Ljr(Q) L 2 (D)). 
Then Eq. (2.1) admits a unique mild solution ( q,r) G CV([0, T]; L 2 (fl; L 2 (D))) xL|(D x 
(0, T); L 2 (D)). hi particular, for any t G [0, T], 


snp E||g(t)|| L 2 (D) + E / \[r(s)\\ 2 L 2 ( D )ds <C< E / ||/(s, •, 0, 0)\\ L2{D) ds + E||g r || L2(D) f, 

ie[o,T] Jo ''Jo J 

where C depends only on A, T and L. Furthermore, if qr G Ljr (Q; Hq(D)), then 

(q, r) e {L 2 ¥ (V x (0, T); H\D ) D Hq(D)) D T|(fi; C([0, T]-L 2 {D))) x L§(D x (0, T); Hq(D)), 

and the following estimate holds 

E( sup ||g(i)||ii ( £»)) + E [ \\<i( s )\\h*(d) + lk(s)llH 0 i ( D) ds 
te[o,T] Jo 0 


<C<E / ||/(s, 0, 0)|| L 2 (D )ds + E||g T || ff i (D ) 


The next lemma provides a standard but useful estimate on the solution to stochastic 
differential equations (SDEs, for short). We list it for ready references. 
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Lemma 2.1. Suppose that x(-) solve the following SDE: 


dx(t ) = (A(t)x(t ) + f(t))dt + (yB(t)x(t ) + g(t))dW ( t ), t G [0, T), 

x(0) = aj 0 , 


where A, B are bounded, /(•) G L 1 ((0, T); R n )) and g(-) G L 2 ((0, T); R n )). 
E sup |x(f)| p < C 

o <t<T 

where C is a constant depending only on A, B. 


EM P + E / | m\dt) +E / \g(t)\ 2 dt 


Then 


3 Approximating BSPDE by BSDEs 


Let {(—be the sequence of eigenvalues and eigenfunctions of A, where {4>i}°L l 
constitutes an orthonormal basis of L 2 (D). is also an orthogonal basis of Hq(D). 

Take the subspace S n = spanjdi, • • • , f> n }, n = 1, 2, • • •. Denote by P n the projection of 
L 2 (D ) onto S n and define A n by A n = P n A\s n . 

The semidiscrete problem corresponding to Eq. (2.1) is to find a process pair ( q n ,r n ) G 
L|( 12; C([0, T]; S n )) x x (0, T); S n ) solving the following BSDE: 

f dq n (t) = ( -A n q n (t) + P n f(t,>,q n (t),r n (t)))dt + r n (t)dW(t), t G [0,T], 

1 q n {T) = P n q T . 


In the next theorem, we prove that (q n , r n ) is convergent to (q, r ) and obtain the rate of 
convergence with respect to the space variable. 

Theorem 3.1. Suppose that (Al) holds and qr G 1/^(0; Hq (D)). Let ( q,r ), (q n ,r n ) be 
solutions to Eq. (2.1) and (3.1), respectively. Then the following estimate holds 


E slip || qit) - q n (t )||| 2(D) +E / || q(t) - q n (t) f H u D) + \\r(t) - r n (t)\\ 2 L2{D) dt 
te[o,T] Jo 0 


< 


C 


(3.2) 


A 


n+1 


M T )\\l^ t ^;H1(D)) + !!/(') 0’ 0)||^2( f2 x(0,T);Hl(D)) 


where C is a constant depending only onT, L and D. 


Proof. We divide the proof into two steps. 

Step 1. Applying Ito’s formula to (q — q n ,q — q n ) l 2 {D)^ f° r an y ^ e [0? T\i we have 

lk(r)-g„(T)||i 2 (D) -||g(t)-^(t)||i 2(D) 

= - 2 f {Aq - A n q n , q -q n ) L2{D) ds + 2 J (f(q, r) - P n f(q n , r n ), q - q n ) L 2 {D) ds 

+ 2 j (q-q n ,r-r n ) L2{D) dW(s) + J \\r - r n \\ 2 L2(D) ds 
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= - 2 


( A(q - q n ),q - q n ) L 2 {D) ds + 2 J ( f(q, r) - P n f(q n , r n ), q - q n ) L 2 (D) ds 

1 r.T 


+ 2 / (q-q n ,r-r n ) L2{D] dW(s) + 

J t J t 

Taking expectation in (3.3), using assumption (Al) 

r>T f*'T' 

(A(q - q n ), q - q n ) L 2 {D) ds 


Ik - r n \\h( D) ds- 

we obtain that 


E ll?(f) ~Qn{t) 1112(D) + E J 


=E\\q(T) - q n (T)\\ 2 L2{D) 
=E\\q{T) - q n {T)f L 2 {D) 


Ik - r n\\mD)< 

rT 


- 2E 


- 2E / ( P, 


L 2 (D)d s 2E 

( f(q, r) ~ P n f(q n , >'„)■ q - Qn) L 2 {D) ds 
( f(s, q, r) - P n f(s , q, r),q- q n ) L 2 {D) ds 
J(s, q, r) ~ Pnf(s, q n , r),q- q n ) L 2 {D) ds 
J(s, q n , r) - P n f(s, q n , r n ), q - q n ) L 2 {D) ds 

n'T 


2E 


2E 


W(I ~ Pn)q\\L 2 (D)\\q - qn\\L 2 (D)ds 


<E||g(T) — qn(T)\\ 2 L 2 ( D ) + 2TE 

J t 

+ 2TE J || (/ — -P n )r ||£2(_D)||g — q n \\L 2 (D)ds 

— Pn)f(s, 0, 0)|| L 2( D )||g — q n \\L 2 (D)ds 
, + \\r - r n \\ 2 L 2 {D) ds 


■2E f \\(I 


-E 


(2 L + l)|k — ?n||| 2 (D) 

^ Ti . 1 JK E IK / - P nW)\\ 2 lU d) + / ||(/ - PnJrWHi.pjds 

>0,0)||^ (D) ds + (4 L + 2)E J ||y - q n \\ 2 tf(p)ds 


J - 

te[0,T] 

•E f \\(I-P n )f{s, 


E j \\r- r n \\ 2 L 2 {D) ds. 
Since for any t G [0, T], 


-2E 


( A(q-q n ),q-q n ) L2(D) ds > 0, 
easily check that 


Jt 

by Gronwall’s inequality, one can 

snp E\\q(t) - q n (t)\\ 2 L 2 {D) 
te[o,T] 

<e (2L+1)T [E||g(T) - q n (T )||| 2(D) + TL max 
L v ’ te[o,T] 


+ TL E ll( 7 - P n)?k)||!2 (D) 


-LE 


|| (/ - P B )r(t) f L 2 {D) dt + E [ T || (/ - P B )/(t, 0, 0) || 2 L 2 {D) dt 

Jo 


(3.4) 


(3.5) 
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Set <?(•) = J2T= i/+•)+■ Then q n (-) = P n q(-) = ELi/+•)+• Hence 


E||9(')-9»OHV)=E|ia-^)9(')fe(D)= E E|^i(-)| 2 < T- E A . E lft 

i=n +1 n+l 


< 


A n+1 




2=1 


< --E||g(.)||^ ( ^ = ^lk(-)ll^ M( D))- 


A n+1 


Similarly, 


(3.6) 


E / ll(7--P»Mt)|||, S(I)) i(i<---E/ ||r(t)||ii, D) <it, 

./n / >n+l 70 

1 


E /_ ll( / -T , n)/(i,0,0)|| L2(D )(it < E I \\f(t,0,0)\\ Hlo(D) dt. 


n+l Jo 


(3.7) 


From (3.4)+3.7) and Theorem 2.1, we have 


< 


sup E\\q(t) ~ q n (t)\\ 2 L 2 (D) 
te[o,T] 
e 3(2L+l)T r 


A 


n+l 


M T )\\l% t (q-,h1(d)) + TL ||g(t)|| 1,2^(£>)) 


< 


C 


A. 


n+l L 


T|l r lli,2( nx ( 0)r ) ; //i(£))) + !!/(•, 0, 0)|| X/ 2(^ x( o, r) ; jF fi( jD) ) 


M T )\\ 2 L^ T (a,H^D)) + ll/(’) 0+)lli|(nx(O,T);fl-i(D)) 


(3.8) 


where C is a constant independently on X n . By the same argument, we obtain that 


E 


< 


C 


A 


n+l L 


I r(t) - r n (t)\\ 2 L 2 (D) dt - 2E / {A(q- q n )(t), (q - q n )(t) ) L i {D) dt 


h( T )\\% T (n-,Hl{D)) + !!/(• + + ) IIl|(Qx(0,T);J71(D)) 


(3.9) 


Step 2. Using (3.3) once more and from similar proceeding as that in Step 1, by the 
Burkholder-Davis-Gundy inequality, we have 

E sup \\q(t) ~ qn(t)\\ 2 L2{D) 
te[o,T] 


=E SJJP] |lk(T)-g„(T)||^ 2 (D ) +2 jf {A(q-q n ),q-q n ) L2{D) ds 
~ 2 [ { f(q, r) - P n f (q n , r„), q - q n ) L2(D) ds 


(3.10) 


- 2 


(q~qn,r -r n ) L2{D) dW(s)- / \\r - r n \\\ 2{D) ds 
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<E||g(T) - q n {T)\\ 2 L 2 {D) + TL max E\\(I - P n )q(t)\\ 2 L 2 {D ) 

+ ||(/ - P n )r(s)\\i, m ds + eJ[ || (/ - P„)f(s,0,0)\\l, m d8 

+ E [ (4L + 2)\\q - q n \\h( D )ds + 2E sup [ (q - q n ,r - r n ) L2(D) dW(s) 
Jt te[o,T]Jt 

<E||g(T) - q n (T)\\ 2 L 2 {D) + TL max E||(/ - P n )q(t)\\ 2 L 2 {D) 

+ LE j '|| (/ - P n )r(s) \\l, (D) dx + E J || (/ - P„)f(s, 0,0) \\l HD) ds 

rT rT -y/2 

+ E / (4L + 2)\\q - q n \\ 2 L 2 {D ) + 6E / (q - q n ,r - r n ) 2 L2 , D) dt 
Jt 1 J o J 

<E||g(T) - q n (T)\\ 2 L 2 {D) + TL max E||(J - P n )q{t)\\ 2 L 2 {D) 

+ LE / ||(/- WS )HW^E [ \\(I-P n )f( S , 0,0) f LHD) i S 


+ (4L + 2)E J ||g - q n \\ 2 L 2 ( D) ds + 18E j \\r - r n \\ 2 L2{D) ds 
+ ^E sup \\q(t) ~ q n (t)\\ 2 L 2 {Dy 

z te[o,T] 

From (3.6)-(3.10) and Theorem 2.1, we obtain 

E^sup^ ||g(t) — q n (t) Hz,2(£>) < ^ ^ \\j( r ^)\\L 2 yrT (n-,H^{D)) + ll/('’0 )IIl2(Qx(o,t) ; ^ 0 1 (d)) 
where C depends only on T, L and D. This completes the proof. 


Remark 3.1. Under the assumptions of Theorem 3.1, applying the results in [22, Section 
4-1], we can prove a better convergence result than that in Theorem 3.1, i.e., 

(q n , r n ) —► (q, r) m L§(Q x (0, T); H\D)) x L|(fi x (0, T); if ^D)). 

However, the convergence speed (given in Theorem 3.1) cannot be improved. 

4 Approximation for BSDE 

In this section, we apply the backward Euler scheme to solve Eq. (3.1). We borrow some 
idea from [10, 29]. 

Since S n = span{0!, • • • ,4> n }, and the solution ( q n , r n ) to Eq. (3.1) is in L[(Q x 
(0, T); S n ) x L|(D x (0, T); S n ), we may take (q n , r n ) to be the following form 

n n 

Qn{') = ^ ] a n,j(')0j) r n(') = ^ ^ , 

1=1 1=1 
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(4.1) 



where a n j £ L|(Q; C([0, T]; R)) and /3 nj - £ L§(fl x (0, T); R). Set 


^n,l (*) ^ 


(fin. l(’)\ 


Ai 

0 

• 0\ 

^n,2(*) 

, fin (') = 

fin, 2(0 

, A — 

0 

to 

• 0 

^n,n(*) J 


\fin,nfi) j 


1° 

0 

• x n ) 


and fnfii (') j fin (')) 


/(/(') 9n(')> r n(')), fil )l 2 (D)\ 
I (/('? 9n(')) r n('))) 02 )l 2 (D) I 


\( /(') Qn('), r n('))i fin )l 2 (D)J 
Then (a n (-), fi n (-)) solves the following BSDE: 


da n (t) = (A a n (t) + f n (t,a n (t), fi n (t)))dt + fi n (t)dW(t), t £ [0,T], 
a n {T) = a nT := ({qr, fix )l 2 (d)i (Qt, 02 )l 2 (d)i ■ ■ • , {qr, fin) L 2 ( D) ) • 


(4.2) 


Before presenting the main theorem, we study the regularity of a n (-) and fi n (-). In the 
sequel, for simplicity, we write d y f n {-, «„(■), fi n (-)), d z f n (-,a n (-),fi n (-)) by d y f n (-), d z f n fi) 
respectively. Similarly, we may use the notations d yy f n (-), d yz f n (-) and d zz f n {-). 

Lemma 4.1. Let (Al) and (A2) hold. Then, for any t,s £ [0, T], 

E\a n (t) - a n (s)| 2 < C( 1 + X 2 n \t - s|)|t - s\, 


where C is a constant depending only on q, T, L and M. 

Proof. Under the assumptions on qr and /, by [10, Theorem 2.6], (a n , fi n ) is in L^ ,2 (R n ). 
Besides, the Malliavin derivative ( Dga n , Dgfi n ) of the solution pair solves the following BSDE: 


0 <8<T 9<t<T 


(4.3) 


Dga n (T) - D e a n (t ) = J^ (AD e a n (s) + d y f n (s)D e a n (s ) + d z f n (s)D e fi n (s )) 

+ J Dgfi n (s)dW(s ), 0 < 9 < t < T; 

l D e a n (t ) = Definft) = 0, 0 <t < 9 <T. 

Moreover, fi n (-) can be represented by D.a n (-). Furthermore, 

sup E sup \D e a n (t)fi + E ( \D e fi n (t)\ 2 dt) < CE\D e a n {T)\\ (4.4) 


J 


where C depends only on q, T and L. 

By above inequality, we have 

E\fi n (t)\ 2 = E\D t a n (t)\ 2 < (E| D t a n (t)fif /q < C(E\D t a n (T)fi) Vq , a.e. t £ [0 ,T\. (4.5) 
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(4.6) 


By (4.2), we also know that 


E|a n (t)| 2 + E/ \/3 n (s)\ 2 ds < C E|a n (T)| 2 + E / \f n (s, 0, 0)\ 2 ds 


Here C is a constant depending only on L and T. By virtue of (4.5) and (4.6), for 0 < s < 
t <T, we can easily have 


E|a n (t) - a n (s)| 2 = E 


(A a n (u) + f n (u,a n (u),/3 n (u))du + / j3 n (u)dW(u)) 


<^4E f A 2 \a n (u)\ 2 du + 16L f (\a n (u)\ 2 + \/3 n (u)\ 2 )du + 8E f f n (u, 0, 0)\ 2 du^j (t — s) 

J s «/ s s 

+ 2E J f3 2 (u)du 

<(AX 2 n (t-s) + 16L(t-s + l) + 8)c(E\a n (T)\ 2 + E J \f n (s, 0, 0)| 2 ds) (t - s) 

+ 2CK\D t a n (T)\ 2 (t-s) 

<C{1 + A 2 (t — s))(t — s ), 


where C depends only on q, L, T and M. This completes the proof. I 

The following lemma is about the regularity of /3 n (-). 

Lemma 4.2. Let (Al) and (A2) hold. Then, for any t,s e [0, T], 

E| f3 n (t) - (3 n (s) | 2 < Ce^-^t - s|(A 2 \t - s\ + 1), (4.7) 

where C is a constant depending only on q, T, L and M. 

We need the following lemma to prove the above result. 

Lemma 4.3. Let (Al) and (A2) hold, mid {T}o<t<T and {4 > }o<t<r which solve the following 
SDEs 

f d&(t) = -V{t) (A + dyf n (t))dt - ^>{t)d z f n (t)dW(t), t e [0, T), 

I *(0) = / n 

and 

d$(t) = (A + dyfnit) + (d y f n (t)) 2 )${t)dt + d z f n (t)$(t)dW(t) t t e [0,T), 

4>(0) = I n , 

respectively. Then, for any p > 2, 

E sup \^(t)\ p <C, E sup \<S>(t)\ p < Ce XnT - 

0 <t<T 0 <t<T 

El sup 4>(s)T(t)| p < C, E|4>(f)T(s)| p < Ce x ^~ s \ 

s<t<T 

E| (<&(£) - 4>(s))T(T)| p < Ce Xnlt ~ sl \t - s\i(\ p \t - s|* + 1), 
where C depends only onp, L and T. 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 
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Proof. First of all, for any Xo € R n , set x(-) = \I/*(-)xo- Then x(-) solves the following SDE: 

f dx(t) — ( - A - d y f n {t)Yx{t)dt - d z f*(t)x(t)dW(t), te [0 ,T), 

\ x(0) = x 0 . 

By Ito’s formula, one can have 


Then, 

therefore 


\x(t)\ p — |x(0)| p 

= f p\x\ p ~ 2 (x, (-A- d y f*x)ds + f p\x\ p ~ 2 (x, -d z f*x)dW{s) 
Jo Jo 

1 

+ -/ X*d z f n d xx \x\ p d z f*xds 

^ Jo 

< [ p\x\ p ~ 2 (x,-d y f n (t)*x) ds + [ p\x\ p ~ 2 (x, -d z f*x) dW(s) 
Jo Jo 

+ 7 :/ X*d z f n d xx \x\ p d z f*xds 
Jo 

<C f \x\ p ds — p f \x\ p ~ 2 (x,—d z ffx) dW(s). 

Jo Jo 


E sup |a:(t)| p < C|a;o| p ; 

0 <t<T 

, , ^, Esup n</<r \x(t)\ p 

E sup |tf(f)| p = sup -<C 

0 <t<T x 0 en n For 


here C depends only on p, L and T. Similarly, one can prove E sup |$(i)| p < Ce 

0 <t<T 

then (4.10) is proved. 


A „T 


and 


Next, we only prove the second inequality of (4.11). The first one can be proved with the 
similar procedure. For any x Q G R n , set x s (t) = < h(t)T(s)a;o. Then x s (t) solves the following 
stochastic differential equation: 

f dx s (t ) = (A + d y f n + ( d y f n ) 2 )x s (t)dt + d z f n x s (t)dW(t), t e [s,T), 

\ x s (s) = Xo- 


By Ito’s formula, one can have 
M^)| p - |x s (s)| p 


= J p\x s (.t )\ p ~ 2 (x s (t), (A dyf n (r) + (d yfn(,T)) 2 )x s (T) ) dr 
+ [ p\x s (r)\ p ~ 2 (x s {T),d z f n x s (r)) dW(r) 
x s (r)*d z f*d xx \x s (T)\ p d z f n x s (r)dT 


1 


<C(A n + l) / \x s (r)\ p dT + C / |x s (r)| p 2 (x s (t), d z f n x s (r)) dW(r). 
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Then, by Gronwall’s inequality, 


E|$(t)^(s)x 0 | p = E|a:,(t)| p < Ce Xn{t ~ s) \x 0 \ p , 
where C depends only on p, L. 

Finally, by Eq. (4.9), one has 
E|($(t) -4>(s))T(T)| p 


=E 


(A + d y f n (r) + (d y f n (r)) 2 )^T)dr^(T) + / d z f n (r)^T)dW(r)^(T) 


<C\ p E( f \*{T)*{T)\dT) P + CE([ |4>(r)T(T)|dry 


GE 


a f n (T)$(T)dW(T)V{T) 


(4.13) 


/ t (*L 

|4>(r)T(T)|dr)" + GE| j d z / n (r)$(r)^(r)$(T) 
:=C(X p n + l)J 1 + CJ 2 . 

For Ji, by (4.11), we have 

./, <E f |$(r)T(T)| p dr( f 1 / E|4>(r)T(T)| p dr(t - s)? 

«/ s J s *J s 

<C(t - s) 1+ 5 = C(t - s) v , 


(4.14) 


where C depends only on p, L and T. 

For J 2 , by (4.11), Holder’s inequality and Burkholder-Davis-Gundy inequality, we also 
can obtain, 


Jo — E 


< E 


a /„(t)4>(t)cHF(t)T(T) 


= E 


2 p\ 1/2 -i /r\ 

a/n(r)$(r)T( S )dlF(r) ) (E|4>( S )T(T)| 2 p) 1/2 


a /n(r)$(r)$(s)^(r)$(s)$(T) 


<C(E( j |$(r)T(s)|W(r)yy /2 (E|4>(s)T(T)| 2p ) 1/2 

1/2 

/ 

<c|(t - s)f J E|4>(r)T(s)| 2p dr} 1/2 (E|4>(s)T(T)| 2p ) 1/2 

1/2 


<C^E 


r-i .!,/■* 
q 


ldrW / |$(r)T(s)| 2p dr) P ] p ) (E|4>(s)T(T)| 2p ) 1/2 


(4.15) 


<G 


S 9 


3 A n(t-s) 


(t-s) } 


<Ge An(t - s) (t-s)2, 


where C depends only on p, L. Combining (4.13)-(4.15), we have (4.12). 
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We are now in a position to prove Lemma 4.2. 

Proof of Lemma 4.2. Let 0 < s < t < T . Since /?„(•) has the representation (3 n (t) = 
D t a n (t), a.s. , a.e. t G [0, T], we have 

Pn(t) - Pn{s) = (. D t a n (t) - D s a n (t )) + ( D s a n {t) - D s a n (s )). (4.16) 


By Eq. (4.3), for any 9 , 6' G [0,t], by Ito’s formula, we can obtain that 

E| D e a n (t) - Dg>a n (t)\ 2 + E J \D e /3 n (s) - D 9 '/3 n (s)\ 2 ds 

=E| D e ot n {T) - Do'Ot n {T) | 2 - 2 (( D e a n (s) - D^a n (s), A(D e a n (s) - D e >a n (s))) 

T ( DgoiniyS') Dg/a n (s), d y f n (s)(Dga n (s) Dgia n (s)) d z fn(s')(Dg[3 n (s') Dg/ /3 n (s )))ds 

<E\Dga n (T) — Dgta n {T)\ 2 + C J \D g a n (s) - Dg^a n (s)\ 2 ds 

+ \ [ | Dg/3 n (s) - Dg,/3 n (s)\ 2 ds, 


where C is independent of n. Here, we apply the fact that A is positive and d y f n (-), d z f n (-) 
are bounded. Therefore, setting 9 = t, 9' = s, by Gronwall’s inequality and (A2) , we deduce 
that 

E| D t a n (t) - D s a n (t )| 2 < CE\D t a n (T) - D s a n (T)\ 2 < C\t - s\, (4.17) 

where C is independent of n. 


For {'h(t)} 0 <t<T and {<!>(£) } 0 <t<r introduced in (4.8) and (4.9), from [28, Theorem 6.14], 
T(t) < h(t) = I n , a.e. t G [0, T\. By virtue of T(-) and $(•), we can list a representation of 
Dga(-). By Ito’s formula, 

d{^{t)Dga(t)) 

= - T(t)(A + d y f n (t))D e a n (t)dt - T( t)d z f n (t)D e a n (t)dW(t ) 

+ T(t)( A + d y f n (t))D e a n (t)dt + ^(t)d z f n (t)D e /3 n (t)dt 
+ V(t)D e p n (t)dW(t) - *(t)d z f n (t)D e Pn(t)dt 
=*(t)(D e p n - d z f n (t)D e a n (t))dW(t). 


Therefore, 


D e a(t )) 

=$(t)T(T)A>a n (T) - Ht) J 4 7(s)(D e /3 n (s) - d z f n (s)D e a n (s))dW(s) 
=Emt)V(T)D 0an (T)\r t ), 
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and 


E| D g a n (t) - D g a n (s)\ 2 

=E\E($(t)*{T)D 0 a n (T)\F t ) - E($(s)tt(T)iV*nCO| F a ) f 
<2E|E(($(t) - <S>(s))-*(T)D e a n (T)\F t ) | 2 
+ 2E\E(^(s)^(T)D e a n (T)\T t ) - E^(s)^(T)D g a n (T)\F s ) 
<2E|($(t) - $(s))y(T)D e a n (T)\ 2 
+ 2E|E($( S )^(T)A ? a n (T)| T t ) - E(^(s)^(T)D e a n (T)\T s ) 

:=2(/i + J 2 ). 

By Holder’s inequality, (A2) and (4.12), we obtain 

h <(E\D e a n (T)\«) 2/q (E(<S>(t) - *(*))¥ (T)|&)^ 
<Ce Xn ^\t-s\(\ 2 n \t-s\ + l), 

where C depends only on q, L , T and M. 

We claim that for any 0 < 9, s < T, 

$(s)y(T)D e a n (T) e M 2 ’ 2 , 


i.e., (T)Dea n (T) has a representation: 


$(s)V(T)Dea n (T) = E (V(T)D e a n (T)) + [ u e>a (t)dW(t ), 

Jo 

and sup e s E| « 0 )S (t )| 2 < C, where C depends only on q, L , T and M. 

If (4.20) is true, then 


h =E E^{s)^{T)D e a n {T)\F t )-E^{s)^{T)D e a n {T)\^ s ) 


=E 


U0 )S (r)dW(r) =E |n 0)S (r)| 2 dr 


<C\t — s|, 


(4.18) 


(4.19) 


(4.20) 


where C depends only on q, L , T and M. Combining the above inequality and (4.16)-(4.19) 
results in (4.7). 


Now we prove (4.20). Indeed, For any 0 < 6, s < t < T, H e ( < h(s)'h(-)) satisfies the 
following SDE: 

' dD g ms)*(t)) = D e (<S>(sMt)){- A - d y f n (t))dt - D e (*(s)*(t))d z f n (t)dW(t ) 

- (^(s)^(t))(d yy f n (t)D e a n (t) + d yz f n (t)D e /3 n (t))dt 
< - Ws)*(t))(d yz f n (t)D ean (t) + d zz f n (t)D 0 p n (t))dW(s ), 9<t<T, 

D e ms)*(0)) = 0 , 

£> 0 ($(s)tf(t)) = 0, 0 <t<d. 
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For any x 0 G R n , set xg tS (-) = D g (ty*(-)$*(s))x 0 and y s (-) = T*(-)<F*(,s)a;o- Then xe, s (-) 
satisfies the following SDE: 

' dx 0 tS (t) = ( - A - d y f*(t))xg ta (t)dt - d z f*(t)xg tS (t)dW(t) 

~ (Dga* n (t)d yy f*(t) + Dg/3*(t)d yz f*(t))y s (t)dt 
< - (Dga^(t)d yz f*(t) + Dg/3*(t)d zz K(t))y s (t)dW(t), 9<t<T, 

xq,s(9) = 0 , 

xe, s {t) = 0, 0 < t < 9 < T. 

For any p G [2, q ), by Lemma 2.1 and (4.11), we have 

E sup \D e ^*(t)^*(s))x 0 \ p 
e<t<r 


<C E / \(D e a* n (t)d yy f*(t) + D e p*(t)d yz f*(t))y a (t)\dt 


y J (4.21) 

/ r T £\ | q _ p 

<C(E sup \D e oi n (t)\ q + e( / |D^ n (i)| 2 di) 2 J (E sup \y s (t)\^) q 

\ 9<t<T ' Jg ' J s<t<T 

P 

<c(supE\D e a n (T)\ q y\x 0 \ p 

<C\x 0 \ p , 

where C depends only on q, L , T and M. Therefore, by the Clark-Ocone-Haussman formula, 

$(s)V(T)D ean (T) 

=E (<f>(s)V(T)Dea n (T)) + J E^D t (^(s)^(T)D e a n (T))y t yW(t) 

=E ms)-*(T)D g a n (T)) 


E(D t (<S>(s)V(T))Dga n (T) + <S>{s)y(T)D t Dga n (T) F t )dW{t) 


:=E (V(T)D e a n (T)) 


ue, s (t)dW (t). 


Fixing p 0 with 2 < p 0 < q/2, by (A2), (4.11) and (4.21), we have 

E|n^(t)| P0 =E E(A($(s)tf (T))A,a n (T) + $(s)tt(T) D t D e a n {T) \7 t ) 

<C(E\D t (^s)^(T))D e a n (T )\ P0 + E\^s)^(T)D t D e a n (T)\ P0 ) 


<C , (E|iL> £ (^>(s)'R(T"))|^h 


E| D 0 a n (T)\ q 


( yy h 

E|<F(s)T(T)|^ 


E| D t D 0 a n (T)\ q 
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where C depends only on q, L, T and M, which proves §(s)'fy(T)Dga n (T) E M 2,P0 C M 2,2 . 

I 


Now, we present the backward Euler method for Eq. (4.2). Suppose a partition n : 0 = 
< ti < ■ ■ ■ < ttf — T of [0, T] with the mesh size 17r| = max^ \t i+ i — tf\. Then we denote 

A i = t i+ 1 - U and A 4 W = W(t i+1 ) - W(ti), for i = 0,1, • • • , TV - 1. 

For simplicity, we assume that A; = 17r| = T, f or each i — 0,1, • • • ,N — 1. Our numerical 
scheme still works for general uniform partition of [0, T] (i.e., there exists a constant K , such 
that K\tt\ < Aj, for any j — 0,1, • • • , N — 1). 

Throughout this paper, we assume that 17r | < 1. 

For the partition n, we introduce the implicit backward Euler method for Eq. (4.2) as 


<(tj+ 1 ) - <(f) =( A <(tj) + fn(tj,<(tj), @n (s)ds | <Ftj ) ) ^ ~ 


rtj+i 


(4.22) 


KMdw(s) 


for any j — 0,1, • • • , N — 1 and a™(t N ) = (gj, (f) l ) L it D y where q J is an approximation of qr- 

Remark 4.1. Multiplying both sides of (4.22) by A jW, and then taking expectation, we 
have 

r t j +1 

E ft(s)ds = E(ci(t j+1 )& j W). 

Jt, 

Furthermore, 


<*n(tj) = A 


-l 

j 


E «(t i+ i)| J).) - f n (tj , al(tj), —E (al(t j+l ) AjW)) 


where 


A; 


/l + XiAj 0 

0 1 + A 2 A j 


° 

0 


\ 0 0 


1 + A n Aj J 


Therefore, this scheme involves the computation of conditional expectations with respect to 
In this respect, the Monte-Carlo method is a popular choice. We refer the reader to the 
related work ([1, 3, 6 ]) for details. 


The following lemma comes from [29, Lemma 5.4], 


Lemma 4.4. Suppose that o* > 0, 6* > 0, c > 0, a t < (1 + cAj)a i+ i + b i+i , for any i = 
0,1, • • • , N — 1. Then 

n 

max at < e cT (a n + > bf). 

0 <i<n 

i= 1 

The following result is the convergence speed for the backward Euler method indicated 
in (4.22). 
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Theorem 4.1. Let (Al) and (A2) hold, and suppose that A^|7t| < 1 is true. Then 

E max | a n (tj) - oT n (tj) | 2 + E [ \(3 n (t) - (3^(t)\ 2 dt 
o<j< N Jo (4.23) 

<CA^(EK(T)-a„(r)| 2 + |x|), 

where C is a constant depending only on Ai, q, T , L and M. 

Proof. We divide the proof into two steps. 

Step 1. Under the condition A^|7 t| < 1, by Lemma 4.2, we have, for any t, s G [0, T ], 

E| Pn(t) - p n (s) I 2 < C\t - s\(X 2 n \t - s\ + 1), (4.24) 


where C is a constant depending only on q , T, L and M. 

For any j — 0,1, • • • , N — 1, taking t = tj in (4.2) and (4.22), one obtains that 
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squaring both sides of (4.25) and then taking expectation, we get 


f-q+i 


E|Aj(a n (tj) -al(tj))\ + E \/3 n (t) - dt 


<(1 + 7Aj/£)E\a n (tj + i) — a^(tj + 1)| + (7 + ej Aj)< E 


rij+i I 2 

A(a n (tj) - a n (t))dt 


4L 2 A? 


H- t-A + L 2 A 2 E| al(tj) - a n (tj )\ 2 + L 2 E / | a n (tj) - a 


9 ' ' J 

rtj +1 


b'+i 


+ Z/E 


+ l 2 e 


E 


I Pn(tj) - A 
ftj+i 


+ L 2 E 
2 


E 


b+i 


An A) ~ Pn(t)dt 




( ^ /3 n (t) - /3 n (tj)dt T t . j | 


<(1 + 7Aj/£)E|a! n (A + i) — a„(A'+i)l" 


fb'+i 


+ (7 + ejA j)< (A 2 A j + L 2 Aj)E / |a n (A) — a n (t)\ 2 dt + 


4L 2 A? 


I 

+ L 2 A 2 E|<(f,) - a n (tj )| 2 + 2L 2 A,E 
rb+i 


9 


fb+i 


| An (A) - An(A| 2rft 


+ L AjE 


mt)-p n (t)\ 2 dt}. 


Therefore, 


1 + AiA j - (j + — ^A 2 ^E\a n (tj) - a*(tj)\“ 



tj +1 


+ H 7 +a-P> |e 

<(1 + 7 Aj / e)E\a n (tj +1 ) — a^(t j+1 )\- 


\f3 n (t) -/3Z(t)\ 2 dt 


ctj + l 


+ (7 + ^/AjO^AnAj + L Aj)E 

4L 2 Af „ r t i +1 
1 2 + 2L 2 A,E 


Set e = 


Ai-7L 2 |vr| 


L 2 


9 

Then 


Ol n (tj ) CTn(t)Adt 

Jt.j 

| An (A) - An(t)| 2 dt 


£ 


7 + )T Aj > 1 — A 1 A ;J > —, 1 + A 1 A ? — ( 7 + 


A, 


L 2 A 2 > 1, 


, >-7 a / , 7 L 2 A j n .. ^ A 1 -7L 2 |vr| 

1 + 7 Aj/£ - 1 + t - vr2 i i> 7 + e/Aj - 7 + 


X 1 -7L 2 \tt\ 


L 2 A, 


(4.26) 


(4.27) 


18 



By (4.27), we can easily obtain that 

i rh+i 

E| a n (tj) - al(tj) | 2 + -E / | p n {t) - Pl{t)\ 2 dt 

Jtj 

7L 2 A ■ 

<(l + ^ _ ^^,\)^\a n (tj + i) — otn(tj+ 1)| 2 + (7 + 


\i — 7L 2 \n\ s f4L 2 Af 


(A 2 A, + L 2 A,-)E 


| oi n {tj) - a n (t)\ 2 dt + 2L 2 AjE / \/3 n (tj) - (3 n (t)\ 2 dt 


Choose 17r| sufficiently small such that Ai — 7L 2 \n\ > Ai/2. By Lemmas 4.1-4.2, we obtain 
that 

i r t i +1 

E\a n (tj) - al(tj )| 2 + -E / |0„(t) - /^(t)| 2 dt 

Jtj 

<(1 + Ci Aj)E|a n (tj +1 ) — a^(tj + i) | 2 

+ + ^A? + 3L 2 A?)(1 + AM) (4.28) 


Ai l 9 


<(1 + CiA^/eI^^+i) - <(ij+i)| 2 + [ \p n (t) - Pl(t)\ 2 dt 


+ C 2 (A n + L +A n |7r|)A J -, 

where Ci and C 2 depend on Ai, q, L , T and M. Therefore, by Lemma 4.4, 

0<j< a iV-2 “ a n( f j)\ 2 + ^E J \Pn(t) ~ 

<e ClT - <(tiv_i)| 2 + ^E / |/3 n (f) -/3£(f)| 2 df 

V 2 -Cat-1 

AT—2 s 

+ C(L 2 + A 2 + Afj7r|) A 2 j 

i=o A 

<e ClT ^E|a n (t 7 v_i) - <(Uv-i)| 2 + ^E \p n (t) - P”(t)\ 2 dt + C(A 2 + A)Vl)M 
Using (4.25) once more, we have 

1 f tN 

E|a n (tAT_i) - a*(t N - i)| 2 + -E / |/3„(i) - Plit)\ 2 dt 


<C(X 2 n + X 4 n \n\)[E\al(T) -a n (T)\ 2 + \n\ 


(4.29) 


(4.30) 


Combining (4.29) and (4.30), we obtain that 


f 1 f tj+1 

max E< a n (tj) - a*(tj) | 2 + -E / |0 n (t) - AT(t)| 2 dt 

o<j<JV-i [ 2 J tj 

<C (xl + AJ|jt|) (e|<(T): - a„(T)| 2 + |tt|) . 
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Now, by (4.28), we have 


N -2 


v _, i 1 /'b+i 

X E|o„fe)-<fe)| 2 + ^E / |/?„(«) -ffl(«)| 2 * 


3=0 

N—2 


'U 


< 22 \ (1 + CiAj)E\a n (t j+ i) — a*(tj + 1)| 2 + C(\ 2 n + A^|7r|)A 

3=0 


Therefore, 

E [ \(3 n (t)-(3:(t)\' z dt = E 


Jo 

N—2 


T , pt N -i r*jv 

7I7+M2 _ nr ( I _j_ 

0 JtN -1 


) \Pn(t) ~ Pn(t)\ 2 dt 


< 


22 2CiAjE(a n (tj + i) — OL^{tj + 1)) 2 + (2 + 2C*iA^v—2)E|cr n (fjv— i) — a n(^iv-i)|" 


i=o 


JV-2 


ftjV 


- 2E|a n (t„) - <(*„) | 2 + C(A, 2 + AJ|7t|) £ A 2 + E / |A,(f) - ®f)| 2 A 

j=0 

<C(A 2 + AIM) (*K(T) - a „(r)l 2 + M 

Step 2. Similar to (4.26), for any j — 0,1, • • • , N — 1, we have 
l A tKfe) - a£(*i))| = |E(Aj(a„(tj) - a£fo)) | ^)| 

fb+i 


<E^|a n (t J+ i) — a^(tj+i)| + 


A(a n (t i ) - a n (t))dt 


+ L 


rfj+i 


A? 4“ l Q; n(^i) a n(^i)l 4" |®n(^j) C^n(t)| 


+ 


1 / /‘ tj+1 

|A»(*i) - Pn(t) | + w-E( / /9£(s) - /3 n (s)ds| Ji. 


+ 


1 

A” 


rb+i 

E( / /3„(s) 


sA 


dt 


Tu 


<E( |a n (tj +1 ) — «n(^'+i)l + 


rtj+i 


A(a n (tj) - a n (t))dt 


+ L A j | a n (tj ) a n ( tj ) | + L 


+ 


E 


rtj+i 


1 / f tj+1 

I Pn{tj) ~ Pn(t) | + ^-E( j PI{S) - /3 n (s)ds|^ 


tj T |c^n(^j) ^n(^)| 


1 

A~ 


b+i 


An (-^) tj n (tj j ds J~ i^ 


dt 


Ft 


(4.31) 
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Hence 


< 


I a n(tj) a n{tj) I 
1 


1 -LA 
+ L 


j 

tj +1 


E( \ain(tj + l) ~ QLn(tj+ 1)| + 


rtj +1 


A () a n (t))dt 


A ,■ 

1 

< 

\1 — L|7t| 

< • • • 


(\A - tj + |a„(tj) - a„(i)| + |A,(fy) -/3„(i)l 

e( /’’ +i 


1 — L|7T 

1 N 2 


E | -I - 6 j 


+ 


Alfa) - 


dt 


Tt 


a n(tj+ 2 ) - «n(^i+ 2 )l|-F tj ) + 


1 — Z/| 7r| 


J+l 


■^) + 7 — 


1 — A 17T 


< 


< 


1 — L|7T 

1 \ AT 

1 — L|tt| 


N-i 

\ N ~i / 1 \ k 

) E (l«n(^iv) - <(ijv)||^tJ + E ( y _ JT~\ ) E ( 6 j+fc-i|*Fii) 

fc=l I 71 "! 


iV-1 


E(|a„M — a£(£jv)| .Ft,-) + 


1 — L|7r| 


"e(^ 6 , 


k=0 


Tt 


<C' 3 {E(|a n (i JV )-<(i JV )||7i J .) +E( J2 b ^)}r 

t k =0 ' 

where C 3 is a constant depending only on T, A; and 

/■h+i rd+i / 

bj =| / A(a n (tj) - a n (t))dt| + A / \V t ~ t i> + K(tj) “ «n(*)l + |£„(tj) - /3 n (*)| 

V t <1 V t A 


. 1 
— E 
1 A; 


h+i 


/^n( s ) _ /5„(s)ds| J 7 ' tj ) | + | —E 


d+i 


Ai(s) - 




sF, I idt. 


Therefore, 


r AT —1 2 ^ 

E 0 ™ Kfe) - <fe)| 2 < C\ E|a n (tjv) - al(t N )\ 2 + e( E &*) (4.32) 

^ V k =0 


For E( we have the following estimate: 


N 1 2 /OF 1 /i /■h+i 


k =0 


e Em E 


3 =0 


A(a n (tj) - a n (t))dt 


rtj +1 


+ L / (y/t - tj + |an(tj) - «n(t)| + |/3 n (t i ) - /3 n (t)|)(Zt 


E (/ ,+1/3 "( S ) -^( S ) d5 l^)| + | E (^ + Vn(s) -^(^1^)1)) 
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N-l 


< 6 N ((A^A j + L 2 Aj)E I ^ | a n {tj) - a n (t)\ 2 dt 

3=0 

r t j + 1 \ r T 

2L 2 A;E / I p n (tj) - p n (t)\ 2 dt + L 2 A 3 ) + 6E / I /3ft) - f3ft)\ 2 dt 


' 2 a 3 

/ l/^nV^V — T -L/ • 

It i ' JO 


N-l 


<C(A 2 + A 4 |tt|) A 2 + 6E / |&(t) - A!(t)| 2 dt 

J =0 ' - y ° 


<C(A 2 + A 4 |tt|) ( E|<(T) - a n (T)| 2 + |tt| ), 
which, together with (4.32), yields that 


E o<i</vI Q!n - C ( X " + Ankl)^K(r) - a n (T) I 2 + (4.33) 


From the assumption, A 2 |7 t| < 1 holds. The above inequality, together with (4.31) and 
(4.33), yields (4.23). This completes the proof. I 

Finally, set 

n n 

= E <<(')*• <(•) = XX‘0*- 

Z = 1 2—1 

From Theorems 3.1 and 4.1, we obtain the following error estimate for our Galerkin numerical 
scheme. 


Theorem 4.2. Let (Al) and (A2) hold, and suppose that A 2 |7 t| < 1 is true. Then 


E o<jxv 11^') ~ &j)\\hiD) + E / IK*) - r n(*)ll! \D) dt 


<C\ 2 n EK(T)-« n (T)| 2 + |7r| + 


C 


A 


n+1 


IkCOU 




+ 


Lg(nx(0,T);Hi(D)) 


Here C depends only on q, T, L, M and D, and olfT) = ( q^, 4>i )l2(d)> where qf is an 
approximation of qr ■ 


Remark 4.2. Generally speaking, one of the main difficulties in constructing a numerical 
scheme for BSDEs is to guarantee the regularity of the second part of the solution. When 
the terminal condition is a function of some forward diffusion, Ma and Zhang [18] obtained 
the L 2 -regularity of the solution's second component, which is the key point of Euler method; 
when terminal condition has no special form, Hu et al. [10] also obtained the L 2 -regularity 
under suitable conditions in terms of Malliavin calculus. In this paper, following some idea 
from [10], we impose assumptions (Al) and (A 2 ). 

In the aforementioned numerical scheme, we suppose that f in Eq. (1.1) is a non-random 
function. As a matter of fact, we can deal with the random case under suitable assumptions 
on f. The reader can refer to [10, Theorems 2.3 and 2.6] for more details. 
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Remark 4.3. In [7], a Galerkin algorithm for the following forward parabolic stochastic 
equation is considered, 


I dU(t ) = (AU(t) + f(U(t)))dt + g(U{t))dW{t), m [0, T\, 
{ U{ 0) = U 0 . 


Under suitable assumption of f, g and U 0 , the discretisation error is bounded by 

C (—^ -1- A^ r+1 / 2 ] + i)|7r| 2r V 

' A n +1 ' 

where 2 r is a positive integer, [x] is the integer part of the real number x and C is a constant 
depending only on U 0 , f, g and T. To some extent our result is consistent with that of the 
forward equations. 
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